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Abstract 


When studying structural vibrations resulting from a concentrated 
source, many structures may be modelled as a finite beam excited by a point 
source. This research explores the theoretical limit on cancelling the 
resulting beam vibrations by utilizing another point source as an active 
controller. Three different types of excitation are considered, harmonic, 
random, and transient. In each case, a cost function is defined and 
minimized for numerous parameter variations. 

For the case of harmonic excitation, the cost function is obtained by 
integrating the mean squared displacement over a region of the beam in 
which control is desired. A controller is then found to minimize this cost 
function in the control interval. The control interval and controller location 
are continuously varied for several frequencies of excitation. The results 
show that control over the entire beam length is possible only when the 
excitation frequency is near a resonant frequency of the beam, but control 
over a subregion may be obtained even between resonant frequencies at the 
cost of increasing the vibration outside of the control region. 

For random excitation, the cost function is realized by integrating the 
expected value of the displacement squared over the interval of the beam in 
which control is desired. This is shown to yield the identical cost function as 
obtained by integrating the cost function for harmonic excitation over all 
excitation frequencies. As a result, it is always possible to reduce the cost 
function for random excitation whether controlling the entire beam or just a 
subregion, without ever increasing the vibration outside the region in which 
control is desired. 
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The last type of excitation considered is a single, transient pulse. A cost 
function representative of the beam vibration is obtained by integrating the 
transient displacement squared over a region of the beam and over all time. 
The form of the controller is chosen a priori as either one or two delayed 
pulses. Delays constrain the controller to be causal. The best possible 
control is then examined while varying the region of control and the 
controller location. It is found that control is always possible using either one 
or two control pulses. The two pulse controller gives better performance than 
a single pulse controller, but the effort to find the optimal delay time for each 
additional controller pulse increases as the square of the number of control 
pulses. 
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Chapter 1 


INTRODUCTION 

With our increasing involvement in the area of aerospace technologies, 
the need to control large, flexible aerospace structures such as aircraft, space 
stations, radio telescopes, or solar energy panels grows more important. 
Unwanted vibrations in these structures may result in poor instrument 
performance, noisy interiors, or even structural failure. 

Traditional methods for controlling unwanted vibrations in a structure 
used passive isolators or stiffened the structure such that it became more 
difficult to excite. Passive isolators are useful only over a small frequency 
band because they are tuned for a single frequency. Isolators also give poor 
performance at low frequencies. Adding stiffness to the structure results in 
increased weight. The design of aerospace vehicles is extremely sensitive to 

weight. 

A novel approach to control, made possible by advances in the digital 
computer, is active control. This method utilizes additional sources of 
disturbance to cancel out the undesired disturbance. Suppressing vibration 
in resonating structures by active control techniques is receiving extensive 
study 1- ®. The plausibility of actively controlling a structure modelled as a 
linear system has been confirmed 7 8 . Very often, the structure to be controlled 
may be modelled as a beam. It may be desired to control the flow of energy 
through the beam or the vibration of the beam itself depending on the 
application. In an infinite or semi-infinite domain, the flow of energy in the 
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beam is the quantity to be controlled. For a finite domain, controlling the 
vibration of the beam itself becomes the problem. This thesis is concerned 
with controlling the vibration of a finite beam. 

If the beam is connecting some mechanical device to a supporting 
structure, vibrations propagate through the beam into the structure. This is 
a situation in which it is suitable to control the flow of energy through the 
beam in an effort to isolate it from the structure without concern for the 
vibration of the beam itself. It has been demonstrated that this flexural 
power flow in beams modelled as infinite or semi-infinite can be attenuated 

9 12 

with a small number of actuators ' 

For finite systems, control is desired over some spatial region. Balas 
describes a method for applying active control to distributed structures and 
demonstrates the method with an example for a simply-supported beam. The 
equation of motion for an undamped beam in bending is solved by the method 
of eigenfunctions. This information is incorporated in a linear control 
system, in which information from a sensor is used to estimate the modal 
amplitudes and velocities. The steady-state controller produces control 
proportional to the modal amplitudes by minimizing a performance function 
via a single actuator. The first three modes of the system after control are 
found to be positively damped for the case of a transient excitation. 

Meirovitch and Baruh 14 also control a simply-supported beam using a concept 
referred to as Independent Modal-Space Control ( IMSC). In this method, the 
modes of the beam are controlled independently, so that one is effectively 
controlling a set of independent second-order systems in parallel. This allows 
for a relatively simple design and quick implementation for real time control. 
The drawback to IMSC is that an actuator is necessary for each mode to be 
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controlled, and in some cases, the number of these can be quite large. The 
large number of actuators makes the control system more expensive and may 
cause design conflicts arising from positioning the actuators. Meirovitch and 
Silverberg 15 demonstrate an application of IMSC by modeling an aircraft 
wing as a cantilever beam and controlling the flutter modes excited by a 
constant velocity airflow. 

All of the studies described above in the area of beam control use some 
scheme to observe the motion of the beam which allows the actuator to 
implement the appropriate control. The sensor provides feedback to the 
actuator so that the vibrating beam is controlled without regard to the input 
disturbance and the actuator continues to damp out the vibrations in the 
beam after the source of disturbance has been discontinued. This thesis is 
concerned with the best theoretical control possible, assuming complete 
knowledge of the beam’s mechanics a priori via the beam equation, and does 
not require any scheme to estimate the state of the system. Since the 
eigenvalues and eigenvectors of the beam are assumed to be known, the only 
additional knowledge that is needed for establishing optimal control of the 
beam is any restrictions on the controller and the input excitation. 

Three different cases of input excitation from a point source on the beam 
are considered, and the resulting vibration is controlled by a second point 
source linearly related to the excitation source. It is this linear relation, 
hereafter referred to as the controller, that is sought such as to give the best 
control possible. 

First the excitation is assumed to be steady-state sinusoidal. A measure 
of control is defined by integrating the mean square displacement over some 
length of the beam. The controller that best minimizes this quantity is then 
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found and its performance demonstrated. The best controller is acausal but 
because the excitation is sinusoidal and therefore completely predictable 
under steady-state conditions, causality is not a constraint for this controller. 
It is referred to as the unconstrained harmonic controller. 

For the second case the beam excitation is assumed to be a stationary, 
white-noise random process. The cost function minimized for this case is 
obtained by integrating the expected value of the displacement squared over 
some length of the beam. The controller that best minimizes this quantity is 
then found and its performance demonstrated. This controller is referred to 

as the unconstrained random controller. 

The last type of excitation controlled is not a steady-state excitation as 
in the case of the previous controllers sought, but a single, transient impulse. 
The form of the controller is chosen a priori as either one or two delayed, 
independent impulses. The delay time and weight of the control pulses are 
chosen such as to minimize the squared displacement response integrated 
over time and some beam length. 

When the best controllers have been found for the previously mentioned 
types of excitation, the results are compared for different configurations. The 
controller location and the interval of control are varied as well the frequency 
in the harmonic case. 
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Chapter 2 


SOLUTION TO BEAM EQUATION 



In this chapter, the equations governing the vibration of a beam are 
developed. The following chapters will use this theory when adding a control 
force to minimize the vibration in the beam. Since the beam equation is 
linear, only one point force need be used in developing the equations of 
motion. Additional point forces are handled by superposition of solutions. 

Given the simply supported, Bernoulli-Euler 16 beam in Figure 1, the 
governing equation for small, transverse vibrations is 


dfyx.t) ( 3*y(x,t) 

dt 2 dx 4 


f(x,t) 


( 1 ) 


with 0 £ x <> l and 0 £ t. The bending stiffness is El, m the linear mass density , 

and f(x,t) the distributed force loading. It is convenient to nondimensionalize 
equation (1) by utilizing the following dimensionless variables for space and 

time: 
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( 1 . 1 ) 



Now equation (1) becomes 


»y(x,r ) + ifyixj) =f{xt) 

dt 2 dx 4 


( 1 . 2 ) 


where 0 ^ x £ 1 , 0 £ t and f(x,t) is a nondimensional loading force given by 


l* 

f{x,t)=Yi f(x ’° 


(1.3) 


For the point force shown in Figure 1, and throughout this thesis 
f(x,t) =f,(t)b(x s ). The simply supported beam is subject to the following 

boundary conditions: 

y(0,f) = y(U) = o 
= o 

dx 2 dx 2 

Equation (2) fixes the ends of the beam while equation (3) ensures that the 
ends of the beam are free to rotate by requiring zero applied moment. 

The homogeneous equation of (1.2) can be solved to find the natural 
vibration frequencies and associated beam modes or eigenfunctions : 

<»„ = (w 7t) w — 1 * 2, 3 t ... 

«J> B (x) = sinnTtx n — 1,2,3,... 

The particular or "forced" solution can be obtained using the method of 
eigenfunctions. Assume an expansion of the displacement in the form 


( 2 ) 

(3) 


(4) 

(5) 


y(r,f)= I ty n {x)y n {t) 

n - 1 


( 6 ) 
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where the generalized coordinate y„(t ) for the nth beam mode is to be 

determined. Substituting equation (6) into (1) and using the orthogonality of 
the beam modes results in the following second-order uncoupled equations 
after adding a modal viscous damping term : 

d\(t) dy n {t) 2 

— tj— + 2 ( 0 n £, —— + (O n y„(t) = 2h(xjf s (t ) ( ■ 1 2 > 

dt dt 

Adding the damping at this point may seem unconventional, but it is a 
common step in the structural field 16 . Putting a damping term in equation ( 1 ) 
would only serve to complicate finding the eigenvalues for the beam. A small 
damping term is included here to prevent an infinite response of the beam at 
resonance. In practice the actual modal damping value ^ is almost always 
determined experimentally 18 . 

Equation (12) can be solved to give 



y„(0 = 2*M f'dx h n (x)f s (t — x) 

Jo 

(13) 

where 

K(0 = ~ e su\(x> d n t u£t) 

toj.n 

(14) 

and 

i-£ 

(15) 


The function h„(t ) is the impulse response of the nth beam mode and u,(t) is 

the Heaviside function. Now with equations (6) and (13) the general solution 
for the steady-state beam displacement can be written as a convolution 
between the exciting force and the impulse response of the beam: 

y(x,t) = f dxh 3 (x,x)f s (t -t) (17) 

Jo 
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where 


( 18 ) 


h s {x,i) = 2 I h n (Wn(x,Wn(x) 

n - I 

The impulse response function h,{x,T) is by definition the response of the 

beam at position x due to an impulsive force at position x, at time x = 0. 

For the case of harmonic excitation, let the force be represented by 

/,(0 = *(//"} 

where /, is in general complex and ft denotes the real part of a complex 

quantity. In harmonic analysis the system is assumed to be in its 
steady-state. In the steady-state, all quantities are oscillating at the same 
frequency «, so all have representations similar to the exciting force in 
equation (19). For simplicity, only the complex amplitudes (indicated by a 
hat) will be utilized in the analysis, and it is understood that the physical 
quantities are found by multiplying by e 1 '* and taking the real part. 

Substitute equation (19) into (17) to give 

y(x,f) = J dx/i,(x,x)9t{/fC y<0< ' _T) } 

Now it can be seen that 

y(x,(0) = H s (x, ©)/, (20) 

where //,(«) is the Fourier transform 19 of the impulse response of the beam. 

The Fourier transform pair is defined by 

H(( 0 ) = j~dth{t)e - ia * < 2U) 

)t(f) = — f d<x>H(<a)e JU * (21.2) 

Fourier transforming equation (18) and (14) gives 
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where 


1 


(23) 


H " ^ (D^-(0 2 + 270X0„C. 

Equation (22) gives the transfer function between the displacement of 
any point on the beam and the exciting force. Figure 2 shows this transfer 
function for the displacement at x = 1/6 due to an exciting force at x, = 1/6 for 
the first 3 modes. The damping was taken to be 1%. The response of the 
modes is well separated and very peaked since the damping is small. The 
phase shift of n radians below and above the resonance frequencies is 
expected. Note the phase change of Jt between the first and second modes. 
This phase change marks the frequency at which the displacement of the 
beam due to the response of the second mode dominates the displacement due 

to the first mode. 

The closed form solution to the beam equation found in this chapter will 
be used throughout the remaining analysis in this thesis. It is worthwhile to 
point out, though, that the remaining analysis is not only applicable to a 
simply supported finite beam. The eigenfunctions and eigenvalues of more 
complicated systems can be found numerically by computer programs such as 
NASTRAN, and utilized in the same manner as those of the beam found in 

this chapter. 
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Figure 2a. Magnitude of the transfer function between exciting force 
and beam displacement at x = 1/6, x, — 1/6, ^, = .01 




Figure 2b. Phase of the transfer function between exciting force and 
beam displacement at x = 1/6, x, = 1/6, = .01 
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Chapter 3 


CONTROLLING HARMONIC EXCITATION 



■ 
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Figure 3. Controller configuration for harmonic excitation 


In this chapter the force exciting the beam is assumed to be harmonic 
with an rms amplitude of unity. A second force is applied to control the 
steady-state vibrations of the beam. The arrangement of the exciting and 
controlling forces is shown in figure 3. The primary source of excitation is 
represented by /(f) and the secondary force, used for controlling the beam, by 
/.(f). The linear relation between the exciting and controlling forces defines 
the controller. In the time domain this linear relation is written 

/ c ( 0 = £^/i ( x)/ ( f-T) ( 24 ) 

where h(t) is the impulse response of the controller. 

For harmonic excitation, the convolution in equation (24) becomes 

/,=//(©)/, ( 25 ) 

where //(<o) is the Fourier transform 19 of the controller impulse response h(t). 
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The displacement of the beam due to a harmonic point source was found 
in Chapter 2. Invoking superposition, the displacement due to both the 
exciting force and the controlling force is 

y(x, to) = H,{x, (o)f,+H c (x,(o)f c {2t 


where 


//,(*, to) = 2 £ $ n (XsWSx) H n(<*>) 

n * 1 

H c (x,(o) = 2 £ ty n (x c )ty n (x)H „(«>) 

n - 1 


( 27 ) 


and 


H n ( ®) = 


I 

CO 2 - G ) 2 + IjdHOnCm 


( 28 ) 


Using equation (25), equation (26) can be rewritten as 

y(x,(s>) = {H,(x, CO) + tf (to)H c {x, to)}/, (2 ' 

Combining equation (27) and (29) gives the harmonic equivalent of equation 


(6): 


Six, a o)= X <j>„(^)y n (w) 

n ■ 1 


( 29 . 1 ) 


where 


y„(o>) =/,{*„(*,) w„(®) 


( 29 . 2 ) 


Note that equation (29) reveals that a controller can be found such that 
the displacement is equal to zero at any single point x on the beam. This 
controller is given by 


H( CD) = - 


H c (x,ti)) 


( 30 ) 


For harmonic excitation, the vibration of the beam will be represented 
by the mean square displacement. Integrating the mean square 
displacement over some length of the beam results in a cost function to be 

mi nimized. Define 
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( 31 ) 


r*i 

<D= dxy\x,t) 

where the bar denotes a time average. The controller that minimizes this 
functional is the controller that reduces the mean square displacement, 
integrated over the specified length of the beam, to the smallest value 
possible with the given beam configuration. This point is stressed because, in 
some cases, this may not be the quantity that is desired to be minimized, and 
if not, the appropriate quantity should be defined. 

Since squaring is not a linear operation, the real part of the complex 
displacement must be taken before squaring: 

y\x, t) = ft 2 {| y(x, to) | (32) 

where 0(jc, to) is the phase angle of y(x, to). Now take the real part and square 


it to give 


y 2 (x,t) =| y (x, to) | 2 cos 2 [tor + 0(.t, to)] 

which can be split into a constant and an oscillating part by the trigonometric 
double angle formula, resulting in 


y\x,t) = l - 1 ;(*,to) | 2 (1 + cos[2tor + 20(jr, to)]} (34) 

Taking the time average of equation (34) and substituting it into 
equation (31) gives 



(35) 


Equation (29) can now be substituted into equation (35), giving the functional 
to be minimized in terms of the sought controller function //(to): 
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<t>(to) = 


I 


dx 


t 


//,(*, to) +//(to)// c (r,to) 



( 35 . 1 ) 


Notice that | f, | 2 12 is the rms amplitude of the sinusoidal excitation which has 


been assumed to be unity so equation (35.1) can be written 


O(to) = 



H s (x, a i) + H (c o)H c (x, to) | 2 


Equation (36) can be expanded in the form 

$(<o) = B,((o) + B 2 *(to)// (co) + B 2 (to)//*(to) + B 3 (to) | H(<a) | 2 


( 36 ) 


( 37 ) 


where 


B,(to) = f dx\H,(xM? 

Jx\ 


C * 1 * 

£ 2 (to) = dxH 1 (x,<£>)H c (x,<n) 

B,(Oi) = ^dx\H c (x,(o)\ z (38) 

Equation (37) cannot be minimized by taking the derivative with respect 
to //(to) and setting it equal to zero because the derivative is not defined for 
complex conjugation. However, the derivative with respect to the real and 
imaginary parts of //(to) is defined. Let 

//(to)=B(to)+./7(to) ( 39 > 

where fl(to) and /(to) are the real and imaginary parts of H (to) respectively. 

Substituting equation (39) into (37) results in 

<D(to) = B,(to)+B 2 (to) (fl(to) + //(to)} +fl 2 (t o) {R((o)-jI (to)} +fl 3 (to) {/? 2 (to) + / 2 (to)} 

( 39 . 5 ) 


Now take the partial derivative with respect to fl(to) and /(to): 
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^ = B m 2 +B 2 + 2B 3 R(G» 
~ = j{B' 2 -B 2 ) + 2B,I((0) 


Set these derivatives equal to zero and solve for /?(©) and /(to). Putting the 
results back into equation (39) yields the optimal controller: 


H( co) = - 


5 3 (<o) 


Using this in equation (37) gives the minimum value of d>: 


<b(co) = B,(o>) - 


I g 2 (o» 1 

B 3 «o) 


(41.1) 


Now use equation (38) to rewrite (41): 


H(i o) = 


dx H,(x, (o)H*(x , (0) 


dx | H c (x,(o) 


It may provide some insight to look at some special cases of equation 
(42). It can be expanded over the modes by using equation (27): 

dx( l> n (jc)<t> m (jr) 

W(<0 ) = -_ — (43) 

1 I H„ ( ( 0 )// * ( co)4» n (x c )<|> m (:t f ) dxty„(x)$ m (x) 

n m 1 m • 1 * x \ 

If d> is to be minimiz ed over the length of the beam, x, =0 and x 2 = 1 . Because 
the eigenfunctions are orthogonal over this interval, equation (43) reduces to 


H{ CD) = 


- 1 \H n (<i>)\ 2 Mx,mXc) 

n - 1 

£ i //„<( o) i 2 <t> 2 (x r ) 
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The frequency response of the controller in equation (43.1) is real, which 
implies that the impulse response of the controller is even 19 . This can be seen 
by writing the inverse Fourier transform 

= (o)e im (43.2) 

271 J— M 


and its complex conjugate, 




( 43 . 3 ) 


Since by equation (43.1), //(to) = //*(»), it can now be seen from equation (43.3) 


that hit) = h(-t). The impact of the even impulse function can be seen by 
looking back to equation (24). To find the controlling force f c (t), future values 
of the excitation force f s {t) are required since h(—t) ^ 0. For this reason the 
controller is referred to as acausal. If the future values of the excitation 


cannot be predicted, an acausal controller cannot be used in a real-time 
application. Obviously the future is predictable for harmonic excitation and 
causality is not a constraint 20 . 

Now that the controller that minimizes d> over the entire length of the 

beam has been found, it is interesting to compare it with the controller that 
minimizes <X> at a single point. This controller can be found by letting the 
interval of minimization approach zero, i.e. | x 2 — r, |— * 0 where the point of 
interest is between r, and jc 2 . In the limit, the integrands in equation (42) can 
be treated as constant and taken outside of the integrations. This results in 


tf(co) = 


-//,(*, (D)//;(.t, to) 

\HXxM \ 2 


( 44 ) 


which simplifies to 
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( 45 ) 


H( < 0 ) = - 


H c (x,<a) 


Equation (45) gives the optimal controller for minimiz ing the mean square 
displacement at a single point on the beam. This agrees with the result of 
equation (30) and from equation (29) it can be seen that the mean square 
displacement at this point is zero. 

Another configuration that can be evaluated analytically is that for 
which the exciting and controlling force are collocated, i.e. 
x s = x c =>H,{x, co) = H c (x, co). It can be seen from equation (42) that H (co) = - 1 
which implies that h(t) = -5(r). Using the definition of h(t) given in equation 

(24), 

m =-j~dxmut - 1 ) = -m m 


which results in perfect control as expected when substituted back into 
equation (36). 

The best controller for harmonic excitation has been found and the 
limiting cases of control checked. Now a digital computer will be used to 
obtain limits on control for more general arrangements. Calculating the 
amplitude of $(x,(o) in equation (29.1) gives the envelope of the beam 
displacement. This amplitude is plotted versus the non-dimensional length x 
of the beam, before and after control, for several excitation frequencies in 
figures 4-7. A table of modal coefficients accompanies each figure. The tables 
reveal the effect of the controller on the magnitude and phase of the 
individual modes. Although all results are presented without dimension, it is 
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necessary to state that the phase will be expressed in radians throughout this 
thesis. Only the first three modes are considered and the coefficient of 

damping for all modes is taken to be 1%. 

In both figures 4 and 5, the driving frequency is at resonance of the 
second mode. In figure 4, <t> is being minimized globaUy or over the entire 
length of the beam, 0 S r S 1 . while in figure 5, it is being minimized locally or 
over a subinterval, 3/4 S x £ 1 . Because the displacement of the beam is 
dominated by the response of a single mode, minimizing the cost function <b 
locally results in very nearly the same global performance as when 
minimi zing globally. 

The contributions of the individual modes when minimizing globally are 
given in table 1. Before control, the second mode clearly dominates the 

displacement. Notice the phase of the individual modes before control. These 

are given relative to the exciting source. The third mode has a phase lag of 
= 0 because it is being excited below resonance; the second mode has a phase 
lag of - n/2 because it is excited near resonance; the first mode has a phase 
lag of - n because it is excited above resonance. After control, the phases are 
unchanged. The reason for this can be revealed by examination of equation 
(29.2). When minimi zing over the length of the beam, H((£>) is real (see 
eqn. 43.1), so it cannot affect the phase of the nth modal coefficient, y„(c »)• 
Assuming a reference phase of zero for the excitation source, the phase of the 

nth modal coefficient is 
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controlled 


Figure 4. Envelope of beam displacement, x s = 1/6, x c = 3/4, v, = 0, 

X 2 = 1, (0 = 0)2 



Before Control 

After Control 

Mode n 

1 ?»(<*>) 1 


1 v„((o) | 

<V„( (0) 

1 

4.839E-04 

-3.136E+00 

1.076E-03 

-3.136E+00 

2 

3.929E-02 

-1.571E+00 

1.749E-05 

-1.571E+00 

3 I 

2.233E-04 

-1.107E-02 

3.600E-04 

-1.107E-02 


Table 1. Modal coefficients before and after control for figure 4. 


_ i [2c«xd„C] 

Z y„((D) = Z//„(to) = -tan' 1 ] — -\ (46.1) 

The modal phases are fixed before a controller is ever found! This explains 
why the phases are unchanged after control in table 1. The amplitude of the 
second mode, which is at resonance, is greatly reduced. When minimizing 
locally in figure 5, the amplitude of the second mode was not as dramatically 
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reduced. This allowed slightly better control within the control interval. 
When the beam was excited at resonant frequencies of the other modes, the 
beam displayed very similar behavior as that seen in figures 4 and 5. 



Figure 5. Envelope of beam displacement, x, - 1/6, x c - 3/4, x, 3/4 

X 2 = 1 , 01 = COj 



Before Control 

After C 

Control 

Mode n 

1 y n m 1 

®) 

l?„(©)l 

^y„( ©) 

1 

2 

3 

4.839E-04 

3.929E-02 

2.233E-04 

-3.136E+00 

-1.571E+00 

-1.107E-02 

1.076E-03 

2.572E-04 

3.600E-04 

1 /V _ J 

-3.132E+00 

-3.105E+00 

-8.591E-03 

i c a K 


Table 2. Modal coefficients before and after control for figure 5. 


Figures 6 and 7 are the same as figures 4 and 5 described above except 
for the driving frequency. For these cases, the driving frequency is between 
the first and second natural frequencies of the beam. Minimizing O over the 
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length of the beam in figure 6 shows virtually no control at all. To 
understand why this is so, use equation (29.1) to expand equation (35) over 
the modes, resulting in 

<D((0) = \ £ £ > n (oi)jC(<o) f dx $„(*)<!>«(•*) ( 4 6-2) 

Z n - 1 m m 1 J jtj 

which simplifies to 

<I>(co) = ^ I I y„(<0) | 2 (46.3) 

4 nm\ 

when calculating <I> over the beam length. This equation reveals that the 

modal amplitudes | _y„(co) | must be independently controlled for global 
reduction of <I>. The table of modal coefficients for figure 6 implies that 
excellent control can be obtained if the coefficients of the first and second 
modal amplitudes can be independently reduced, but this is not possible 
because the controller //(co) has only one degree of freedom when minimizing 
globally (see eqn. 43.1). 

Between resonances, two modes are contributing on the same order to 
the total response of the beam. Therefore, controlling one mode excites the 
other even more, and this spillover limits the amount of control that is 
possible when minimizing globally. When minimizing locally, the controller 
uses spillover to achieve better control as seen in figure 7. Notice that the 
beam displacement outside the subinterval of control is much larger than 
before the control was applied. Table 4 gives the contributions of the 
individual modes to the displacement envelope in figure 7. The amplitudes of 
all three modes are increased significantly, but the phase of the second mode 
was changed so as to cancel the contributions of the first and third mode in 
the interval of minimization which created a larger displacement envelope 
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outside the interval of minimization. It will be seen in figure 14 that moving 
the location of the controller allows better control both inside and outside of 
the minimi zation interval. 



Figure 6. Envelope of beam displacement, x, — 1/6, x c - 3/4, x, 0, 
x 2 = 1, 01 = 2.5(0, 



Before Control 

After C 

Control 

Mode n 

1 £.(©) 1 


1 y„(co) 1 

<v„(co) 

1 

2 

3 

1.382E-03 

1.289E-03 

1.942E-04 

-3.132E+OOm 

-2.050E-02 

-6.019E-03 

1.120E-03 

1.488E-03 

1.758E-04 

-3.132E+00 

-2.050E-02 

-6.019E-03 


Table 3. Modal coefficients before and after control for figure 6. 


Now that the beam displacement envelope has been examined for 
several cases and the amount of reduction over the length of the beam seen, it 
is interesting to study the reduction of d>, which corresponds to the 
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Figure 7. Envelope of beam displacement, x, = 1/6, x c = 3/4, v, = 3/4, 
x 2 = 1, (o = 2.5co, 



Before Control 

After Control 

Mode n 

i y„(o) i 

^v„(< o) 

,iy„«o)| 

<v„( CO) 

1 

1.382E-03 

-3.132E+00 

8.614E-03 

-3.061E+00 

2 

1.289E-03 

-2.050E-02 

4.226E-03 

-3.052E+00 

3 

1.942E-04 

-6.019E-03 

7.019E-04 

5.494E-02 


Table 4. Modal coefficients before and after control for figure 7. 


performance of the controller. Without control, it can be seen from equation 
(37) that d> = 5,(co). After applying control, <D is given by equation (41.1). 
These two quantities are compared in figures 8-13. 

Figures 8 and 9 give the performance of the controller versus frequency 
for minimizing over the beam length, and over a subinterval of the beam 
length, respectively. Figures 4-7 become single data points on these graphs. 
The largest reduction of d> is at the resonance frequencies in both figure 8 and 
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9. Notice that there is a frequency between the modes at which no control is 
possible. The exact frequency at which this occurs is a function of a, and -V 
At this frequency there is a balance of the modal responses in the region of 
control, so spillover prevents any reduction of the cost function. The location 
of Figure 9 shows that a greater reduction in <D is possible when minimizing 

over a subinterval of the beam length. 



Figures 10-13 illustrate the importance of location of the controller if 
that freedom is available. All of these plots show the perfect control when the 
controller is collocated with the source. The beam is being driven at the 
resonant frequency of the second mode in figures 10 and 11. Note that 
insignificant global control is possible in figure 10 when the controller is 
located at a node of the mode that is dominating the displacement, but figure 
11 shows that local control is always possible at a resonance, no matter 
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where the controller is Iocs. ted. 

In figures 12 and 13 the driving frequency is between the first and 
second natural frequencies. If one is trying to control globally, figure 12 
shows that it is best to be as close as possible to the source. When 
minimizing locally, a more complex relationship results, as shown in figure 
13. Interestingly enough, the best control is achieved (for separated sources), 
when the controller is located at the center. At this position the controller 
cannot excite the second mode and as a result avoids spillover which allows 
greater freedom in altering the first and third modes, and results in much 
better control. This is the location of the controller in figure 14, in which all 
parameters are identical to those in figure 7, except for the controller 
location. It can be noted from table 5 that the second mode is not excited by 
the controller, because the amplitude and phase are unchanged, after control 
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uncontrolled 

controlled 


X c 

Figure 10. Effect of controller position x c on reduction of global cost 
function d>, x, = 1/6, x, = 0, x 2 = 1, © = ©2 



Figure 11. Effect of controller position x c on reduction of local cost 
function d>, x, = 1/6, x, = 3/4, x 2 = 1 , (O = o >2 





The amplitudes of the first and third modes are increased and the phase of 
the first mode is shifted by n, so that the optimal control is achieved by a 
destructive interference between the modes in the region of minimization. 



Figure 12. Effect of controller position x c on reduction of global cost 
function d>, x, = 1/6, x x = 0, x 2 = 1 , to = 2.5(0, 
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Figure 13. Effect of controller position x c on reduction of local cost 
function <I>, x s = 1/6, r, = 3/4, x 2 = 1 , to = 2.5(0, 

Figures 7 and 14 illustrate control of the beam over an interval. In 
figure 7 control was possible at the cost of a considerable increase in 
displacement outside of the interval of minimization. In figure 14, much 
better control was shown to be possible with smaller displacements outside 
the interval of minimization with the proper location of the controller. For 
these same two locations of the controller, it is interesting to examine the 
effect of a continuous change in the interval of minimization. Comparing the 
global cost function to that obtained when minimizing locally (over a 
subinterval), gives an indication of the efficiency of the controller. To do this, 
fix the right side of the interval at x 2 = 1 and let 0 < x t < 1. The ideal controller 
is then found from equation (41) and d> from equation (37). 
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Figure 14. Envelope of beam displacement before and after local 

minimization with optimal controller placement, a, = 1/6, 
x c = 1/2, A, = 3/4, x 2 = 1 , CO = 2.5(0, 



Before Control 

After Control 

Mode n 


<V„((D) 

1 y n (a>) 1 


1 

1.382E-03 

-3.132E+00 

1.441E-03 

-2.339E-02 

2 

1.289E-03 

-2.050E-02 

1.289E-03 

-2.050E-02 

3 

1.942E-04 

-6.019E-03 

3.925E-04 

-1.451E-02 


Table 5. Modal coefficients before and after control for figure 14. 


Figure 15 shows almost no increase in the global cost function when the 
interval of minimization is changed. This indicates that minimizing locally is 
the same as minimizing globally. When any single mode dominates the total 
response of the beam, then minimizing the response of that mode results in 
both a global and local minimum. The result is much different if the exciting 
source is off resonance, as can be seen in figure 16. In this case, the global 
cost function 0 increases dramatically to allow greater minimization in the 
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(0 = (Dj, x, = 1/6, x c = 3/4, x 2 = 1 

interval of concern. This is a result of spillover being utilized to reduce O 
locally. If the controller is located optimally, this increase in the global cost 
function O, due to spillover, can be minimized as shown in figure 17. It can 
be seen that with the proper location of the controller, tremendous reduction 
in <D can be achieved locally without a significant increase in the global value 

of O. 

In summary, the best control of harmonic excitation over the length of 
the beam is achieved when the displacement is dominated by one mode. Over 
the length of the beam, the modes must be independently minimized and a 
single controller can only independently minimize a single mode. If control is 
desired only over a subregion, it is possible to control all modes over this 
subregion if spillover outside of the region of control is allowed. In general, 
the best location for the controller is as close as possible to the exciting 
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CO = 2.5©,, JC, = 1/6, x c = 3/4, x 2 = 1 



Figure 17. Effect of interval of minimization on both the global and 

local value of cost function <t> before and after control with 
x c optimal, ©=2.5©,, x s = 1/6, x c = 1/2, x 2 = l 
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source. Results similar to these presented have been seen in experimental 
control of a semi-infinite beam 9 , which would seem to indicate that 
comparable mechanisms are present for other types of boundary conditions. 
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Chapter 4 


CONTROLLING RANDOM EXCITATION 



In this chapter the force exciting the beam is assumed to be a random 
process with a constant power spectrum and a variance of unity. This type of 
random process is also referred to as white noise because of its uniform 
frequency content. A second force is applied to control the stationary 
vibrations of the beam. The arrangement of the exciting and controlling 
forces is shown in figure 18. The primary source of excitation is represented 
by f a (t) and the secondary force, used for controlling the beam, by f c (t ). The 
linear relation between the exciting and controlling forces defines the 
controller. In the time domain this linear relation is written 

f c (t) = j~dxh {x)f s (t - x) (47) 

where h(t ) is the impulse response of the controller. 


33 



The displacement of the beam due to a single exciting force is given by 
equation (17). By superposition, the total displacement is 

y(*»0= I dT,A # (x,T 1 )/j(f-T l )+ f dX x h c (x,X x )f c (t-X x ) (48) 

Jo Jo 

Using equation (47), this can be written 

y (*,<) = l dx x h s (x,x x )f t (t-x x )+ dx x h c (x,x x ) )j dxh(x)f J (t-x l -x ) (49) 


For random excitation the vibration of the beam will be represented by 
the expected value of the displacement squared. Integrating this quantity 
over some length of the beam results in a cost functional to be minimized 
with respect to the impulse response function of the controller. Define 



(50) 


as this cost functional where, o^(r) is the variance of the displacement which 


is identical to the expected value of the displacement squared. To find T in 
terms of h(t), the impulse response of the controller, first square equation (49) 

y\x,t)= f dx x dx 2 h 3 (x,x x )h 1 (x,x 7 )f s (t -x x )f s (t -x^ 

Jo 

+ 2 J ^dxh (x) dx x dx 2 h 1 (x, x x )h e (x , tj/,(f - x x )f s (t -x 2 -x) 

+ jjxdx'h(x)h(x') J o dx x dx 7 h c (x,x x )h c (x,x 2 )f s (t-x x -x)f s (t-x 2 -x') (51) 
and then calculate the expected value, 
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°5c*) = j dx i d \ R f(^ - h)K(x> \)h,(x, x,) 

Jo 

+ 2 Jj/tA(t)J^ dX x dX 2 R f (X x — Tj — T)/j,(Jt, T,)/t r (jC, T 2 ) 

+ ^T 1 rfT 2 /? / .(T,-t 2 +t-T')/l r (jr,T I )/l f (A,X 2 ) (51.1) 

where ^(x) is the autocorrelation of the random process For white noise 

with unity variance, R/ t (x) = 5(x). Using this in equation (51.1) and integrating 
with respect to Xj yields 

oJ(jr)=r dx x h 2 s {x, x,) 

Jo 

+2 J <#xA W j[ dx^Cr.x^/t^x.-x) 

+J dxdx' A(x)/i(x / ) dx,/»^.r,T l )/t c (jr,x 1 +x-x') (51.2) 

The cost functional T is found by integrating equation (51.2) over some 
length of the beam: 

¥ = &,(()) + 2 J dxh(x)b 2 (x) + j dxd't h(x)b i (x-x')h(x r ) (53) 

where c m C x ^ 

Wjo dx \ J d xh s (x,x x )h 3 (x,x x -x) 

^( T ) E Jo <f**,( x » T i)M*.*|-X) 

00 r x 2 

i/x, dxh c (x,x x )h c (x,x x -x) (54) 

0 J jr, 


/ 
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The Fourier transform of &,(x), b 2 ( x), and & 3 (t) in equation (54) results in £,(©), 

B 2 ((o), and fi 3 (co), respectively. These transforms have previously been 
determined in equation (40). Now use the definition of the inverse Fourier 
transform defined by equation (21.2) to find T in the frequency domain: 


I d(oB,((o) 

2icJ— 

P ** 1 C °° y(o> | + 

+2 J </t— J dco^coj/^co^^co^e 

+J d (0, doij d (Oj// (co,)B 3 (o) 2 )//(co3)c 

Using the identity 19 


;(m, + 


(55) 




(55.1) 


equation (55) can be rewritten as 

T = — f ~dco (5,(0)) + 2£ 2 *(co)/f (co) +5 3 (w) I H (co) | 2 } (55.2) 

2tw-~ 

From equation (53), notice that B’ 2 (<o)H (co) must be a real quantity, so 


equation (55.2) can be rewritten as 

T = — f **d<» {B ,(<») + Bl((o)H (co) + 5 2 (co)//*(co) + fl 3 (<o) | H (co) | 2 } (56) 

2jcJ— 

From equation (37) it can be seen that the integrand of equation (56) is <t>(co), 
so it can be written 

T = — f d co<D(co) ^57) 

2k 

It is apparent from equation (36) that 4>(co) is an even function, so equation 
(57) can be written 
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T = - f </(0<D(to) 

K Jo 


(58) 


With equation (36), equation (58) can be expanded into 

'F= f dx-( d(£>\H^x,(x))+H((i))H c (x,()i)\ 2 (59) 

Jx | K Jo 

which, when compared to equation (50), gives 

afc) = - f d(d | H 3 (x, (0) +H(w)H c (x, co) | 2 (60) 

K Jo 

To minimize 'F in equation (58), O(co) must be independently minimized 

at each value of co, because the function O(co) is, by definition, positive for all 
to. This independent minimization has already been accomplished in chapter 
3. The unconstrained controller that gives the min imum vibration for 
harmonic excitation is also the best unconstrained controller for random 
excitation! After some thought, this fact is not so surprising. Without the 
constraint of causality, the impulse response of the controller is defined for 
<>t <,<*>. This has the effect of letting the controller "know" the input 
excitation before it excites the beam. Thus the random input is not random 
as seen by the controller, but deterministic. The infinite limits in equation 
(47) illustrate this mathematically. In this light, it makes sense that the 
optimum controller is the same for both random and periodic excitation. 

Now that the unconstrained controller for random excitation is known, 
the performance of the controller will be examined. It is interesting to note 
that the special cases outlined in chapter 3 also hold for the case of random 
excitation; perfect control is achieved when the exciting source and controller 
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are collocated, and when minimizing the response at a single point on the 
beam. Because these cases result in O(co) a 0 for all values of (o, equation (58) 
verifies that the same holds for controlling random excitation. 

For more general cases, the square root of equation (60) gives o v (.v), the 

deviation of the beam displacement. It is plotted against x to show a 
statistical beam displacement envelope in figures 19 and 20. In all plots, only 
the first three modes are considered and the coefficient of damping for all 
modes is taken as 1%. In figures 19 and 20, the vibration is being minimized 
from 0 ^ x <, 1 and from 3/4 < * < 1 , respectively. The first mode dominates the 
uncontrolled displacement envelope in figure 19, but this is hardly surprising 
since the relative dominance of the first mode was seen in figure 2. These 
results are similar to those realized before for harmonic excitation. 

Compared to the global minimum, better control was achieved in the region 
of minimization, at the cost of less control outside of the region of 
minimization. 

The vibration of the beam over the interval of minimization is 
represented quantitatively by 'F in equation (58). In figures 21 and 22, T is 
plotted versus the position x c of the controller. Perfect control is achieved for 
collocation of source and controller. When minimizing over the length of the 
beam in figure 21, the reduction in T lessens as the controller position 
coincides with nodes of the second and third modes. This is a result of the 
requirement that each mode be minimized independently for global reduction 
of T, and a mode cannot be controlled if it cannot be excited. The 
requirement to min imize the modes independently can be verified by 
substituting equation (46.3) in equation (58) to expand T over the modes: 
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Figure 19. Deviation of beam displacement, 0/x) before and after 
global control, x, = 1/6, x c = 3/4, x x = 0, x 2 = 1 



Figure 20. Deviation of beam displacement, o v (x ) before and after 
local control, x. = 1/6, x c = 3/4, x, = 3/4, x 2 = 1 




v p= 7- 1 f ^wiy n (co)i 2 

47tn -1 Jo 


( 61 ) 


Even though the controller cannot control a mode that it cannot excite 
because of actuator location, it can still control the remaining modes, so 
figure 21 shows a mini mum of 15 dB in the reduction of the cost function 'F. 
This is quite different than what was observed in the case of harmonic 
excitation, where some controller location always existed for which no control 
was possible. 



Figure 21. Effect of controller position x c on global cost function T , 

x s = 1/6, x t = 0, x 2 = 1 

When minimisin g over a subinterval of the beam, figure 22 illustrates 
that locating the controller near a node can result in improved control. 
Figure 24 illustrates this improved control by plotting the deviation of the 
displacement after locating the controller at x e - .49. Because of the modal 
coupling in a subinterval, the first and third modes are utilized to locally 
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The effect of the interval of minimization on the reduction in 4* is seen in 

figure 23. Greater control is possible as the interval of minimization is 
decreased. The global value of T is also reduced when controlling locally for 
all values of the lower limit of minimization. 



Figure 24. Deviation of beam displacement, o y (x ) after optimal 
placement of controller, x, = 1/6, x c = .49, x x = 3/4, .v 2 = 1 

The performance of the unconstrained controller has been examined, 
and it establishes that excellent control is possible over a broad frequency 
bandwidth for random noise excitation using only a single controller. 
Regardless of the interval of minimization, a reduction in the global beam 
vibration is always achieved. This differs from the results of chapter 3 
where, for harmonic excitation, an increase was observed in global beam 
vibration when minimizing locally for certain arrangements of exciting 
frequency and controller location. The best location of the controller for 
global control was shown, in general, to be as near the exciting force as 
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possible, while avoiding node locations of the modes to be controlled because 
each must be controlled independently. When controlling locally, however, 
positioning the controller at a node location can result in improved control 
because of the local interference between the modes. 

The controller was revealed to have the same transfer function as the 
controller used for harmonic excitation. This controller is not constrained by 
causality which may limit its use depending on the application, so a 
controller constrained by causality is also of interest. The discussion of the 
constrained random controller is taken further in appendix A. 
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Chapter 5 


CONTROLLING TRANSIENT EXCITATION 



In this chapter the force exciting the beam is assumed to be an impulse 
function. Two different controllers are used to control the resulting transient 
vibrations. The first controller uses one impulse delayed in time to minimize 
the vibrations, whereas the second controller uses two impulses, both delayed 
in time. The delays constrain the controller to be causal. The equations will 
be formulated for the two pulse controller, allowing the single pulse 

equations to be found by letting the amplitude of the second controlling pulse 
vanish. 

The arrangement of the exciting and controlling forces is shown in figure 
25. The primary source of excitation is represented by f s (t) and the secondary 
force, used for controlling the beam, by/ r (r). For a single pulse excitation and 
two pulse control, these can be written 

(69) 
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( 70 ) 


m =/,. 8 (' - 0 + /c2^ - h) '2 > ^ 0 

The linear relation between the exciting and controlling forces defines 
the controller. In the time domain this linear relation is written as 

f c ( 0 = J dxh(x)f s (t -x) (7 

where h(t) is the impulse response of the controller. Substituting equations 
(69) and (70) into (71) gives 

h{t) = h l S(t-t ] ) + h 2 b(t-t 2 ) ( 


where 


h - — 


h - — 

K -f, 


sO 


( 73 ) 


The displacement of the beam due to a single exciting force is given by 
equation (17) where h£x,t) is defined as the response of the beam at position x 
due to an impulse located at x,. By superposition, the total displacement is 

y (x , t) = J ~dx h,{x , x)f s (t - x) + jjlxh c (x,x)f c (t - x) ( 74 ) 

where h c (x,t) is defined as the response of the beam at position x due to an 


impulse located at x c . Using equations (69) and (70), this can be written, 

y (x, t) = f s0 h s (x, t) +f c iK( x > f ~ + fc 2 ^ci x ^ ~ r z) 


( 75 ) 


Normalize by the amplitude of the exciting impulse by defining 


y s (x,0 


y (*>0 

fs O 


( 76 ) 


which allows equation (75) to be written 

y 8 (x,r) = h s (x , 0 + h x h € {x,t - r,) + h 2 h c (x , t - 1 2 ) 


( 77 ) 
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In the case of a single pulse controller, lu-0 and equation (77) is left with 

two degrees of freedom, h, and r„ Similarly, the two pulse controller has four 
degrees of freedom. The four degrees of freedom will be reduced in the 
following analysis to three by letting the time between the two pulses be fixed 

at some positive constant a such that a = t 2 - t x . 

A measure of the transient vibrations resulting from the initial impulse 

will be formed by squaring the displacement and integrating over all time. 
Integrating this quantity over some length of the beam results in a cost 
function to be minimized with respect to the impulse weights and delays of 

the controller. Define 

A.prJVvW) < 78 


To find A, first square equation (77): 

y\{x, t ) =h]{x,t) + 2h x h s (x, t)h c (x,t - t x ) + 2 h 2 h,(x,t)h c (x,t - t 2 ) 

+ h\h]{x,t - t x ) + 2h x h 2 h c {x,t - t x )h c (x, t-t 2 ) + h\h]{x ,t -t 2 ) (79) 

Substituting equation (79) into (78) gives the cost function to be minimized 
explicitly: 

A = A, + 2h x A 2 (t x ) + 2h 2 A 2 (t 2 ) + h 2 A,(t x , f.) + 2h l h 2 A,{t l , t 2 ) + h 2 A,(t 2 , t 2 ) 
where A x = ['dxj~dth 2 (x,t) 

A 2 (x,)s f 'dx f dth s (x,t)h c (x,t-x ,) 

Jx s •'O 

A x ( t„t 2 ) = f 'dx f dth c (x,t-x x )h c (x,t-x 2 ) 

Jx. JO 


( 80 ) 


( 81 ) 
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It can easily be seen that the definitions in equations (81) are the same as 
those already defined by equations (54) where A t = 6,(0), A 2 ( x,) = b 2 { x t ), and 
^3(^0 T2) — ^3(^2 ~ ^1) = Now equation (80) may be written as 

A = 6,(0) + 2 6,6 2 (f,) + 2 h 2 b 2 (t 2 ) + 6 2 6 3 ( 0) + Ihfob^a) + 6 2 6 3 (0) (82) 


To find the minimum value of the cost function A, calculate the partial 
derivatives with respect to the unknowns h u 6,, and t,. Setting these 
derivatives equal to zero yields stationary points, which are necessary 
conditions for a minimum. The derivatives are 


dA 

3^ =2/1,63(0) + 26 2 6 3 (a) + 26 2 (r,) 

(83) 

3^ = 2h M<*) + 26 2 6 3 (0) + 2b 2 (t 2 ) 

(84) 

dA 

3^ =26,6' 2 (/,) + 2h 2 b' 2 (t 2 ) - 26,6 2 6' 3 (a) 

(85) 


Setting equations (83) and (84) equal to zero yields the linear system 


f> 3 (0) 63(a) 

L^a) 6j(0) 


W = |-6 2 (f,)l 

W 1 - 6 2 (f 2 )J 


which can be solved to give 


_ ^2^1)63(0) 6 2 (f 2 )6 3 (a) 
6 3 2 (0) - 6 3 2 (a) 


6 , = - 


6 2 = - 


62(^63(0) -6 2 Q,)63( a) 
63(0) - 6 3 2 (a) 


( 86 ) 


(87) 


( 88 ) 


After setting equation (85) equal to zero, equations (87) and (88) can be 
used to write a single equation in r, whose roots are stationary points of A. 
However, the equation is transcendental, with infinitely many roots. 
Moreover, it is not clear which root results in the global minimization of the 
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cost function A. An easy way around this problem is to substitute equations 
(87) and (88) into equation (82) to give A as a function of r, and then to graph 
the cost function so that the global minimum may be identified visually. This 

is done in figures 26 and 27. 

In figure 26, the cost function A has been minimized over the entire 
length of the beam, or globally. Both the one and two pulse controller achieve 
their minimum values at nonzero delay times. The minimum for the one 
pulse controller is reached at t x = T,/2 - .318, where T x is the period of the first 
eigenvalue, or eigenfrequency, of the beam, equal to 2 /tc = .637. The delay 
times at which no control is possible for a single pulse controller correspond 
to the instants at which the energy in the first mode is all potential. Any 
attempt to control the beam with a single pulse at this instant only increases 
the energy in the first mode. It is possible to control the higher order modes 
at these instants, but spillover into the dominant first mode (see figure 2) 
prevents any reduction of the cost function. 

To better illustrate this mechanism, consider a single degree of freedom 
spring mass system in an excited state. As the mass reaches the system 
equilibrium position its energy is purely kinetic and an impulse can be used 
to completely stop the subsequent motion. This is in contrast to the time at 
which the main mum displacement is reached and the energy of the system is 
purely potential. At this time, the best controlling pulse would be none at all, 
because striking the mass with an impulse can only serve to increase the 
energy in the system. The situation is more complicated for multiple degrees 
of freedom but the mechanisms are the same. 
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Figure 27. Cost function A vs. delay time t x of initial control impulse 
for local control, x s = 1/6, x c - 3/4, x, = 3/4, x 2 = 1 
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The single degree of freedom example in the last paragraph explains 
why the fundamental period of the oscillating cost functions in figures 26 and 
27 are half the period of the first beam mode, but there are several other 
trends that can also be explained. One of these is the upward trend of the 
local minima for the cost functions. Recall that the plotted cost function is 
integrated over all time (see eqn. 78), but it is an impulse that is being 
controlled, so the major contribution to the integral comes from the time 
interval immediately following the exciting impulse. Thus, to minimize the 
cost function it is, in general, better not to delay control too long, because 
after the exciting impulse has decayed away, there is nothing left to control 
and the oscillating cost functions will have converged to the uncontrolled 

level. 

Again looking at the single pulse control in figure 26, it is seen that the 
global minimum for A occurs at a delayed value for r, rather than a zero 
delay, and the fourth minimum results in better control than the third. It 
seems that waiting later can result in better control of the beam. This is true 
because of the interactions of the beam modes. The best delay at which to 
control is the delay at which the spillover from controlling the first mode 
helps control other modes as well. This phenomenon appears to be a strong 
function of the controller location as seen in figure 28, in which the optimal 
delay time /, for the global cost function is plotted against the controller 
location. The changes in optimal delay times are discrete jumps that occur 
when the controller location is such that the controller spillover couples with 
the beam modes in a way to make an alternate delay time more desirable. 
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Figure 28. Effect of controller position on the optimal delay time for 
the controller, x, = 1/6, r, = 0, x 2 = 1 , a = T t /4 

In contrast to the single pulse controller, the two pulse controller in 
figures 26 and 27 always results in some reduction of the cost function. This 
results from the controller being able to offset the spillover into the first mode 
with the second pulse while controlling higher order modes. The amount of 
control possible is highly dependent on the time a between the two pulses. 

For increasingly small values of a, the cost function ceases to show any 
noticeable changes after values smaller than <x = T./lOO, but better control is 
possible for greater choices of a. The increased complexity of the cost 
function for the two pulse controller will be discussed later. 

It is worthwhile to comment on the similarity of figures 26 and 27 since 
some may be surprised that a global reduction would yield very nearly 
identical values oft , to minimize the cost function A. This is a result of the 
transient nature of the excitation. For transient excitation, no steady 
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opposition to a local disturbance can be developed. For this reason, the global 
and local minimiz ations yield very nearly the same results, as will be shown 
in later plots. 

Now it is of interest to note the controller weights /r, and h 2 that yield the 

control shown in figures 26 and 27. Figures 29 and 31 show these weights 
corresponding to the global control in figure 26, and figures 30 and 32 show 
those for the local control in figure 27. All of the weights oscillate with the 
period of the first beam mode. 

For a = r,/ 4 in figures 29 and 30, the first weight of the two pulse 

controller is essentially the same as that of the single pulse controller; the 
second weight is the same also except that it is delayed by the time lag a 
between the control pulses. The peak values correspond to the optimal delay 
times in figures 26 and 27 and the zero crossings coincide with the delay 
times at which no control is possible. 

The amplitude h x of the single pulse controller at f, = 0 can be largely 

attributed to the effect of the excitation source and controller locations on the 
first mode. To better see this, assume the beam is represented by its first 
mode only. Minimizing the cost function A in equation (78) with respect to 

= 0) results in /i, = sm( 7 tx,)/ sin(7Ut c ). For the locations of source and 
controller in figure 29, /i, =sin(jc/6)/ sin(3re/4) =1/^2 ~ .707. The actual value of 
.66 read from figure 29 is lower because there is more than one mode to 
control and stopping the first mode completely by letting h, = -11^2 would 
result in spillover into the remaining modes, resulting in less than optimum 
control. 


52 



1 



0.6 H 


1 

<D 

5 


* 

2 . 


O 

c 


0.4 H 


0.2 - 


0 - 
- 0.2 - 


- 0.4 - 


- 0.6 - 


- 0.8 - 


-1 - 

0 


Figure 29. 


+ hi , 1 pulse 

ti hi , 2 pulse 

h 2 , 2 pulse 


Weight of controller impulses for one and two pulse global 
controllers, x, = 1/6, x c = 3/4, r, = 0, x 2 — 1» ot = T x i 4 



Figure 30. Weight of controller impulses for one and two pulse local 
controllers, x, = 1/6, v = 3/4, a, = 3/4, x 2 — 1, <x= TJA 




In view of the slight difference between figures 26 and 27, it is no 
surprise that figure 30 is very similar to figure 29. The increased importance 
of the higher order modes in figure 30 accounts for the increased distortion of 
the oscillating values of the weights. 

Now look at the case in which the time a between the controlling pulses 

gets small. The weights become large and most often opposite in sign, as 
shown in figures 31 and 32 for a = 7,/ 100. The "dipole" nature of the pulses 
allows control of higher order modes while limiting spillover in the first mode. 
The peak values in these figures indicate the delays at which the single pulse 
controller could not achieve any control. This two pulse controller does 
achieve some control (less than 3 dB globally) at these peaks as seen in 
figures 26 and 27, but larger control pulse weights are needed. The 
magnitudes of these pulse weights will be associated with the effort of the 
controller and are given in figures 31 and 32. The amplitudes are over ten 
times the amplitudes for more widely spaced controlling pulses that not only 
need much less effort but result in better control as well. 

The best control results near the zero crossings of the weights in figures 
31 and 32. Here the two weights sum to that of the single controller as given 
in figures 29 and 30. This agrees with the results seen in figures 26 and 27, 
where all configurations of one and two pulse controllers gave identical 
performance at the value f, = T,/2. Thus, in return for somewhat increased 
performance over a single controller pulse, the two pulse controller requires 
significantly more effort if a « 1 . 

Figure 26 showed a large change in the optimum delay time r, of the two 
pulse controller when the time a between the control pulses was changed. To 
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Figure 31. 


Weight of controller impulses for two pulse global 
controller, r, = 1/6, x c = 3/4, r, = 0, - 1, a - T,/ 100 



Weight of controller impulses for two pulse local 
controller, x, = 1/6, = 3/4, x, = 3/4, r 2 = 1 , a = r,/ 100 


Figure 32. 




better exhibit this relationship, a contour plot of the reduction in the cost 
function A is shown in figures 33 and 34 with t, and a as independent 
variables. For a = 0, the determinant of the linear system in equation (86) is 
zero because t, = t 2 , so results for the single-pulse controller are plotted along 
the ordinate for this degenerate case. 

There is more to be seen in figures 33 and 34 than that they are complex 
and contain many local minima for A. Both figures contain strong contour 
lines that can be explained with the help of figures 26 and 27. The diagonal 
contours are centered about the two lines, f, + a = T,/2 and t, + a = 37" 1 /2, which 
correspond either to the first or second pulse acting at the minima seen in 
figures 26 and 27 for the single pulse controller. The strong vertical contours 
are located about the lines f, = T t /2 and r, = 3T./2 which again correspond to the 
minima seen in figures 26 and 27. Along these contours, the second pulse has 
virtually no impact at all, indicating that its amplitude must be near zero. 

The global minima in figures 33 and 34 occur near the origin and show a 
reduction in the cost function A of more than 16 dB both globally and locally. 
The minima of A are at (r„ a) = (.06, .12), and (>„ a) = (.07, .1 1) in figures 33 and 
34, respectively. The location of the minimum is caused by an interaction 
between the modes. Note the general spreading of the contour lines as r, and 
a increase. This can be understood by recalling again that the cost function A 
is obtained by integrating over time, so that, in general, the longer the 
controller delays the controlling pulses, the less reduction is possible. 

Recall that the contours in figures 33 and 34 are functions of the 
controller configuration. To study the effect of controller location and 
interval of minimization, it would be necessary to regenerate and search 
contours such as figures 33 and 34, for each individual location of the 
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controller and interval of minimization, so that the values of f, and a that 
result in global reduction of the cost function A would be found. Because this 
is time consuming, a is held fixed in figures 35-39, and the optimum t , sought 
along the line a = One should keep in mind though, that the 
quantitative control is also a function of the time a between control pulses. 

Figures 35 and 36 show how the location of the controller affects the 
quantitative control. Just as in the case of controlling a harmonic exciting 
source in chapter 3, collocation of source and controller results in perfect 
control. It can be seen from figures 35 and 36 that substantial control is 
possible for any location of the controller in the case of global minimization of 
A, but, in general, it is best to move the controller as close as possible to the 
exciting source. Again a very subtle difference is seen between the shape of 
the global and local curves. This is seen throughout the plots in this chapter 
and reinforces the earlier statement that minimizing A locally in the case of 
transient excitation also results in very nearly a global minimum as well. 

Now that the effect of controller location on A has been illustrated, it is 

of interest to see at what effort the control was gained. The effort can be seen 
by examining the weights /i, and h 2 of the control pulses. The weights 
corresponding to the global control in figure 35 are shown in figure 37 and 
the weights corresponding to the local control in figure 36 are shown in figure 
38. When the controller location nears the ends of the beam, the weights get 
large because poor modal coupling makes it increasingly difficult to excite the 
beam modes. The sudden changes in magnitude of the weights is due to a 
change in the optimum delay time r, for the initial controlling pulse. These 
jumps were seen earlier in figure 28. As the controller position is changed, 
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Figure 35. Effect of controller position on global cost function A, 
x s = 1/6, =0, x 2 = 1, a= jT,/4 
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Figure 36. Effect of controller position on local cost function A, 
x s = 1/6, r, = 3/4, x 2 = 1 a = T,/4 





the optimum value for f, will change in discrete jumps. Correspondingly, the 
weight of the controller pulses will also change in discrete jumps. The sign 
change of the single pulse controller results from the change in sign of the 
modal eigenfunctions as the controller location changes. The second mode 
dominates this change in sign, which causes the jump to be near x=\!2 which 
is the location of the sign change for the second modeshape. When either of 
the two weights h x or h 2 equals zero, there is no advantage to using a two 
pulse controller. These regions match the regions in figures 35 and 36 where 
no further reduction in the cost function A was possible with the two pulse 

controller. 

The consequence of changing the controller location x c was seen in 

figures 35 and 36. The location x c has a major impact on the amount of 
control possible. Figure 39 shows that changing the interval of minimization 
does not have any dramatic results on the amount of control possible. 

The effect of changing the controller location, interval of minimization, 
delay time of the controlling pulses, and the time between the control pulses 
have been observed. The performance was measured by the reduction in the 
cost function A. Recall that A is integrated over time and the interval of 
minimization. To better understand physically how the cost function is being 
reduced, define a spatial energy density function such that 
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Figure 38. Effect of controller position on weight of control pulses for 
local control, x s — 1/6, x, = 3/4, x 2 = 1 a = T,/4 




(89) 


A = 1 dxE yW 

where r « 

£ »W S J 0 dty\x,t) (90) 



The function E y (x) shows where on the beam the reduction in the cost 

function A is taking place. It is plotted in figure 40 for both global and local 
control by the one and two pulse controllers. The delay times t { and a in 
these plots are the optimal values as selected from figures 33 and 34. 

The plots for the single pulse controller overlay one another and are 
indistinguishable. The interval of minimization has virtually no effect on the 
control possible by a single pulse. Looking back at figure 39 confirms that, 
regardless of the interval of minimization, the same global control is achieved 
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pulsed, = T,/2), 2 pulsar, = .07, a = .1 1 )), globally, = 0, x 2 = 1 , 

1 pulsar, = T,/2), 2 pulse(f, = .06, a = .12)} 

with a single pulse controller. This is not true of the two pulse controller. It 

gives better control in the interval of minimization at the cost of giving up 

some control outside the interval. 

From figure 40 one gains a feel for the location on the beam at which the 
cost function is being reduced. A more physical demonstration of how the two 
pulse controller achieves better control than the single pulse controller is 
evident in figures 41 and 42. These show time histories of the beam 
displacement y(x,t) at a single position on the beam before and after control. 
The two pulse controller again uses the optimal values for r, and a as selected 
from figures 33 and 34. 
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Figure 41. Time response of beam before and after control at position 
x = 3/4 using a single pulse controller, x, = 1/6, x, = 0, x 2 = 1 , 


uncontrolled 

controlled 


Time response of beam before and after control at position 
* = 3/4 using a two pulse controller, v = 1/6, x, = 3/4, .v, = 1, 





Figure 41 gives the displacement before and after a single pulse 
controller is used. The before and after control plots are identical until the 
controller strikes at delay time The first mode is virtually eliminated from 
the time response after application of the controller. Notice the delay of the 
controller corresponds to r, = T,l 2 which is a time of zero displacement of the 
first order mode. At this instant the energy in the first mode is entirely 
kinetic which, as mentioned in the discussion of figures 26 and 27, marks the 

best time to strike the beam with a controlling impulse since the first mode 
dominates the displacement. 

The two pulse controller in figure 42 strikes sooner than the single pulse 
controller in figure 41 in an attempt to reduce the displacement of the first 
beam mode in its initial half cycle. This is accomplished, but at the cost of 

slightly increased oscillations following the control pulses as compared to the 
single pulse controller. 

In summary, it is always possible to control the transient vibration of 
the beam with both one and two impulses. The two impulse controller results 
in increased performance, as seen in figure 40, but at the cost of increased 
effort in the determination of the delay time and weight of the second pulse. 

To make the weight of the second pulse simple to determine, the time a 
between the pulses can be fired as it was for most figures in this chapter but 
this offer, only a slight improvement over the single pulse controller. Thus, 
for the two pulse controller to offer significant improvement, as was the case 
in figure 40, it is necessary to find the optimal delays for the first and second 
pulse independently. This is a point of difficulty for this method. Roots of a 
transcendental equation yield infinitely many choices for the delays ,, and a 
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that result in stationary points of the cost function A, but there is no clear 
method for quickly determining the values of t , and a that result in the global 

minimum . 
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Chapter 6 


SUMMARY 

A controller to optimally reduce the vibration of a finite beam has been 
found for each of the three cases of excitation considered. For each controller, 
the performance was illustrated for a variety of configurations and the 
important results will now be summarized. 

For harmonic excitation the optimal controller is seen to have two 
degrees of freedom, the amplitude of the controlling force and its phase 
relative to the exciting force. When minimizing the vibration globally, 
however, the only degree of freedom of the controller is the amplitude of the 
controlling force, because the orthogonality of the modes requires that each 
mode be minimized independently. Perfect control is possible at any single 
point on the beam, but when minimizing the vibration over some length of 
the beam, the only way to achieve perfect control is to collocate the excitation 
source and the control source. For separate locations, the best control exists 
when the excitation frequency is near resonance of one of the modes. 

Between the modal resonances, the quantitative control diminishes, until no 
control is possible at a discrete frequency which is a function of excitation 
location, controller location, and the interval of minimization. This was true 
when minimizing both locally and globally, but the frequency band in which 
the control was diminished, was much broader when minimizing globally. 

This bandwidth in which little control is achieved should not be a foremost 
concern though, because it always occurs at frequencies where the response 
of the beam is at a minimum and it is not as important to control the 
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vibration. The location of the controller should generally be as close as 
possible to the excitation source. For global control, avoid locating the 
controller at a node of a mode to be controlled; when controlling locally, 
locating the controller at a node of a contributing mode can be advantageous 
because controller spillover into that mode is eliminated. 

For random excitation, the ideal controller was found to be the same 
controller that best minimizes the vibrations resulting from harmonic 
excitation. The performance of the random controller was then shown to be a 
frequency average of the performance of the harmonic controller, resulting in 
excellent control of random excitation using only a single controller. Perfect 
control is possible at any single point on the beam, but when minimizing the 
vibration over some length of the beam, the only way to achieve perfect 
control is to collocate the excitation source and the control source. The best 
location of the controller for global control was shown, in general, to be as 
near the exciting force as possible, while avoiding node locations of the modes 
to be controlled, because each mode must be controlled independently. When 
contro lling locally, however, positioning the controller at a node location can 
result in improved control because of the coupling between the modes. 
Regardless of the interval of minimi zation, a reduction in the global beam 
vibration is always achieved. More control is gained as the interval of 
minimiz ation is reduced, but not at the cost of increasing the global energy. 

For tr ans ient excitation, two different controllers were found and 
compared. The first controller used a single impulse delayed in time to 
minimize the vibrations, whereas the second controller used two impulses, 
both delayed in time. The delays constrain the controller to be causal. If the 
time between successive pulses was fixed for the two pulse controller, it 
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resulted in virtually no improvement over the single pulse controller. The 
transient nature of the excitation prevents the controller from being able to 
zero the vibrations at a single point on the beam as was possible for 
steady-state excitations. Perfect control is achieved if the excitation force and 
controlling force are collocated. Changing the interval of minimization had 
no effect on the single pulse controller, and only a slight effect on the two 
pulse controller; Both controllers exhibited excellent control over all intervals 
of minimization. The best location for the controller was again as close as 
possible to the excitation source. 

The highlights on controlling beam vibrations resulting from three types 
of excitation have been presented. The study was by no means exhaustive, 

and there are several interesting topics that could be examined further as 
future research. 

In the case of both harmonic and random excitation, the controller was 
found to be have the same acausal transfer function. To implement an 
acausal controller, future values of the excitation source are required. This is 
not a problem in the case of harmonic excitation, since all future values are 
exactly known, but exact future values of a random excitation are generally 
not available. This requires the practical constraint of causality to be 
imposed on the controller used in the case of random excitation. It would be 
interesting to find the transfer function of the optimal controller constrained 
by causality and to compare its performance to the acausal controller in 
chapter 4 to see what is lost in performance as a result of the constraint. 

It would also be useful to extend the theory of controlling a single 
transient pulse given in chapter 5 to control a sequence of random pulses. 

The controller need not be limited to only one or two control pulses. If many 
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pulses are used and the time between the control pulses is fixed, this problem 
is exactly a discretized version of that of finding the optimal controller 
constrained by causality for random excitation mentioned previously. 

Also it would be desirable to implement and verify the controllers found 
in this work digitally in an experimental arrangement. It is not necessary, 
though, to use the controllers found in this work. Adaptive controllers could 
be used to control the vibrations, and their performance compared with the 
theoretical limits presented in this thesis. 
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Appendix A 


THOUGHTS ON CAUSAL CONTROL OF A 
STATIONARY RANDOM PROCESS 

In chapter 4, an optimal controller was found to control white noise 
excitation. The disadvantage of the controller was that it was acausal and 
required future values of the noise excitation in order to be implemented. In 
some applications, a good estimate of future excitation values may exist, but 
for a completely general application, the controller for random excitation 
must be constrained to using only past values of the random noise excitation. 
Unfortunately, the constraint of causality makes the ideal controller 
considerably more difficult to find. The calculus of variations is employed 
here to find an integral equation that the constrained controller must satisfy. 
Both the analytical and numerical solution of this equation are then 
discussed. 

Imposition of the constraint of causality is possible in both the frequency 
and time domains. In the frequency domain, causality causes the real and 
imaginary part of the transfer function to be dependent. If the real part is 
known, the imaginary part can be found by the relation 21 



One way to satisfy the constraint of causality is to let //(co) = R((a) + jl(Oi) in 

equation (56), where /?((o) and /(to) are the real and imaginary part of H(t o), 
respectively. Now use the integral relation in equation (Al) to get the cost 
functional 'F in terms of only one independent function. This results in 
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which is a very complicated integral equation. 

In the time domain, causality is more simply enforced by requiring h(r), 

the impulse response of the controller, to be zero for t < 0. This affects the 
limits in equation (53) such that 

4 , = &,(0) + 2 f dTh(i)b 2 (x)+ ( dzdY h(x)b 3 (x-?)h(t') (A2) 

Jo Jo 


Now the problem is in the time domain rather than the frequency domain. 
Unfortunately, the quantity to minimize is still an integral equation, but not 
as complicated as was the frequency domain equation. 

An impulse response function h(t) is sought that will minimize T, which 

is by definition, always a positive quantity. This is a problem in the calculus 
of variations 22 . If a minimum exists, a necessary condition is that it be at a 
stationary point of ¥. Finding a stationary point is not sufficient however, to 
say that the point represents a minimum or maximum, but fortunately in 
many physically motivated problems, it will be obvious. If only one 
stationary point exists for T, then it must be the minimum since a controller 
could always be chosen to make T larger. To find a stationary point, first 
calculate the variation of T in equation (A2) with respect to the controller 
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m- 

54^ = 2 J dxb 2 (x)hji + J dxdx'{b 3 (x-x')h(x')8 x h+b 3 (x-X’)h(x)S,h} (A3) 

and rearrange to give 

84* = 2 f dxb 2 {x)hjx + f dxdxf {b 3 { x - 1 7 ) + b 3 (xf - 1 )} h (x')S r /i ( A4) 

Jo Jo 

Looking back at equation (54), b 3 (x) is 


r x i 

b 3 (x) = j dx x J dx h c (x,x l )h c (x,x i -x) 


(A5) 


Now use X 2 = x, - 1 to change the variable of integration resulting in 




-M 


dx h c (x,x 2 + x)h c (x,x 2 ) 


(A6) 


which, when compared with equation (A5), clearly reveals that b 3 (x) = b 3 {-x). 


Because b 3 { x) is an even function, equation (A4) can be rewritten as 

54* = 2 J dx^b 2 (x) + j dx'b 3 (x-X')h(x'^b T h (A7) 

A stationary point is defined by 84* = 0. Since 5 is a function that can 

be chosen arbitrarily to prevent the integration from being zero, the 
integrand itself must everywhere be identically zero. This results in the 
Euler, or governing, equation: 

b 2 (x)+ f dx'b 3 (x-x')h(x') = 0 (A8) 

Jo 

This particular integral equation is of a very common form known as the 
Wiener-Hopf integral equation. If the controller were not restricted to t > 0, it 
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could be solved by utilizing the convolution theorem of the Fourier transform 
which gives the same result as already found for the unconstrained random 
controller in chapter 4. 

It may seem that the properties of the one-sided Laplace transform could 
be utilized to give 


H(s) = - 


BAs) 
B 3 (s ) 


(A9) 


but this is not possible since the convolution theorem is defined only for 
functions that are zero for t < 0; b 3 (x) is even, so the Laplace transform cannot 
be used. Of course the integral in equation (A8) could be separated into two 
one-sided integrals but then h(x) is coupled with limits and an algebraic 
solution such as equation (A9) is not possible. 

An analytical method for solving equation (A8), which requires the 
Fourier transform of the kernel b 3 ( x) to be factored into two functions that are 
analytic in the left and right half-planes, respectively, is presented by 
Krein 28 . Factoring the product of two summations over the modes would be a 
very tedious process, but it may be possible to utilize a computer program 
that does symbolic manipulation, such as Macsyma 24 , to aid in this task. 

A straightforward approach to obtain a numerical solution is to 
discretize equation (A8) and solve the resulting simultaneous linear system of 
equations for discrete values of the desired controller function. Discretizing 
equation (A8) yields 

N 

b 2 (kA) + A L b 3 (kA-nA)h(nA) = 0 k =0, 1, ...,N (A10) 

n -0 
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The functions b 2 {i) and b 3 (t) are presented in closed form in appendix B. The 

solution to equation (A10) yields the impulse response of the controller. 
Substituting this into a discrete version of equation (A2) gives the 
approximate performance of the controller: 

N N N 

'¥ = b l (0) + 2A X h(nA)b 2 (nA) + A 2 £ 1 h(nA)b 3 (nA-kA)h(kA) (All) 

rt-0 * -o * -0 

Unfortunately, the solution of the discrete system created questions rather 
than answers. The discrete values of the impulse response h(t) had large 
oscillations at t = 0 and t = N A. The oscillations became more tightly spaced 
and stronger as the stepsize A was decreased. Substituting the result into 
equation (All) indicated a large reduction in the beam vibration. To check 
the reduction of displacement in the time domain, the beam equation was 
discretized and control simulated. After control, the amplitude of the 
displacement increased! If the oscillations were removed from the ends of the 
discrete impulse response, a reduction in displacement was observed. This 
would seem to indicate that the oscillations were not part of the actual 
solution, but if they were eliminated, then other problems occurred As an 
example, minim i z ing the beam vibration globally did not reduce the global 
cost function as much as minimizing locally did. This makes absolutely no 
sense at all. This may be the result of the determinant of the system matrix 
rapidly approaching zero as N increases, although no thorough explanation 
could be found. 

The problem of finding the optimal causal controller has been briefly 
discussed. This is a problem which contains many subtleties, and further 
research is necessary before these subtleties can be fully understood. 
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Appendix B 


RELEVANT INTEGRALS 


The closed form solution of several integrals is needed for plotting the 
solutions found in the text. The first is 

h) = I dx h n (x - t x )h m { x - r 2 ) t 2 > t x £ 0 

JO 


where 


h n(t) = ~~e^'sii \<£> Jn t u s (t ) 


co, 


a, n 


This can be solved to give 

h) = COS COj J t COStoj JzX", 

-cose sin <a dm t 2 Y m 
- sin C0j n r, cos <a dm t 2 Y mn 
+sin(x> d j i sine o d ^ m t 2 Z im } 


where 


and 


co rf , n eo Jm 

x sin ( t0 ^- + ^.m)h ~ (c A& + co^Cm) cos(eo rf n + eo ,,Jr 2 

2{(eo„C, + cn„Cn) 2 + (eo,,. n + eo^.J 2 } 

ft 0 '.- - oyj sin(eo,, n - eo,,, m )f 2 - (eo n ^, + eo,^) cos(eo,, n - eo^ m )e 2 
2{(to n ^, + eo m ^,) + (eOj „ — eoj,,,) 2 } 

Y J3& + + (co^.n ~ eo,,, J cosjw^, - eo,, Jr 2 

" 7 2{(eo n C n + co m C n ) 2 + (co,. n -eo^ 

, (gj. + m X) sin(eo,, „ + eo,,, m )r 2 - (eo^, n + eo^ w )cos(eo^ n + eo,,, Jr 2 
2{(co n ^ + eo m ^,) 2 + (eo,, „ + eo,,^) 2 } 
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(tty „ - (o^J sin(o> rfi „ - (Oj >m )r 2 - (&>„£. + co^ ^Jcostoy,, 
Znm ~~ 2{(co n ^ + (o m ^) 2 + ((o Jn - Mj.J 2 } 

((Oj t „ + (Q I^J sin(co^ „ + COj m )f 2 - ((OnC. + ^mCw) 008 ^" + 
2{((0„£„ + ®«C") 2 + (®J." + 

Another integral needed throughout the text is 

C m (x v x 2 ) = f dxty n ( x)ty m (x) 

Jx t 

where <J> n (jf) = sinnnx n = 1,2,3,... 

The solution is 

sin(n — m )kx 2 — sin(n - m )7tx, sin(n + m )xx 2 — sin(n 4- m )7ir l 


C nm (x v x 2 ) = 


2 K(n -m) 


2jt(n +m) 


= - (r, - x.) - - (sin 2 tj 7EX 2 - sin 2/ncx,) n-m 

2 4nn 


n 
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